Orbital-selective Mott-Hubbard transition in the two-band Hubbard model 
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Recent advances in the field of quantum Monte Carlo simulations for impurity problems allow 
-within dynamical mean field theory- for a more thorough investigation of the two-band Hubbard 
model with narrow/wide band and SU(2)-symmetric Hund's exchange. The nature of this transition 
has been controversial, and we establish that an orbital-selective Mott-Hubbard transition exists. 
Thereby, the wide band still shows metallic behavior after the narrow band became insulating -not 
a pseudogap as for an Ising Hund's exchange. The coexistence of two solutions with metallic wide 
band and insulating or metallic narrow band indicates, in general, first-order transitions. 
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By virtue of dynamical mean field theory (DMFT) m 
, our understanding of the Mott-Hubbard transition 3] 
in the one-band Hubbard model has greatly improved in 
the last years. The bandwidth-controlled Mott-Hubbard 
transition is, at least within DMFT 0,0], of first-order 
at low temperatures (T) and becomes a smooth crossover 
for temperatures above a critical point, which terminates 
the first-order line. A further complication arises exactly 
at zero temperature where two solutions coexist like for 
low Ts. But at T = 0, the insulating solution is always 
higher in energy than the metallic one, i.e., the insulating 
solution is metastable throughout the whole coexistence 
region. The DMFT Mott-Hubbard transition is of second 
order at T = despite the coexistence of two solutions. 

For making contact with experiments, orbital realism 
has to be taken into account, e.g., within the merger of 
local density approximation and DMFT (LDA-^DMFT 
approach |5|). In the case of transition metal oxides, 
typically either the three <2g or the two Cg bands cross 
the Fermi energy. At the very least, these orbitals should 
be included. For degenerate orbitals, the situation seems 
to be clear, at least within DMFT: there is a first-order 
Mott-Hubbard transition 6J. Most transition metal ox- 
ides are, however, non-cubic. Hence, the orbital degen- 
eracy is lifted. Take, for example, the unconventional 
superconductor Sr2Ru04 l7| which has a wide d^y band 
and narrow dxz,zy bands [Sj and which becomes a Mott- 
Hubbard insulator upon substituting Sr by Ca 

For such a situation with wide and narrow bands the 
details of the Mott-Hubbard transitions are so far in- 
conclusive, even within DMFT and even for a simple 
two-band Hubbard model with Coulomb interaction U 
and Hund's exchange J between the two bands: Koga 
et al. employed the so-called exact diagonalization 
(ED) method to solve the DMFT equations and obtain 
two Mott-Hubbard transitions: first the narrow band be- 
comes insulating, then the wide band. This scenario 
has been coined orbital-selective Mott-Hubbard transi- 
tion In contrast, Liebsch |0| uses quantum Monte 
Carlo (QMC) simulations and the iterated perturbation 
theory (IPT) to solve the DMFT equations and finds 



a single first-order Mott-Hubbard transition with simi- 
lar changes in both bands. On the insulating side, the 
wide band has a pseudogap which gradually amplifies to 
a real gap with increasing U. In principle, the QMC 
is more suitable for addressing the Mott-Hubbard tran- 
sition since ED only gives discrete peaks in the spectra, 
making it difficult to unambiguously identify a gap. How- 
ever, the QMC simulations are restricted to relatively 
high temperatures and there is a sign-problem 13] if the 
Hund's exchange coupling is taken into account in full, 
i.e., not only the Ising but also the spin-flip component. 
Since the same limitations as in |l2l| also apply to all 
previous LDA-HDMFT(QMC) calculations |3|, there is 
an urgent need to clarify whether and how the details of 
the Mott-Hubbard transition are affected. Another im- 
portant aspect is whether two solutions coexist. Liebsch 
finds two coexisting solutions at a single transition, while 
Koga et al. do not address this question. If there was 
a first-order transition two consecutive transitions might 
even be bridged into a single one. 

In this paper, we study this transition by employing 
the most recent advances in the field of QMC simu- 
lations for DMFT. The recently introduced projective 
QMC (PQMC) method Q enables us to address T^O. 
Furthermore, the new Hubbard-Stratonovich decoupling 
of 15] allows for the calculation with the full SU(2)- 
symmetric Hund's exchange at a still-manageable sign- 
problem, in particular in combination with PQMC. 

Model. Starting point is the two-band Hubbard model 
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Here, c\^^ and c^^^ are creation and annihilation opera- 
tors for electrons on site i within orbital m and with spin 
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a. The first line describes the kinetic energy for which we 
employ the semi-elliptic non-interacting density of states 
^°(^) = /8 V(W^m/2)^ - £^ (orbital-dependent hop- 
ping amplitudes 1^ on a Bethe lattice). For the following 
calculations, we use different widths for the two bands: 
W\ = 4 for the wide and W2 = 2 for the narrow band as 
in 12] . Note that there is no hopping/hybridization 
between the two bands. The second line describes the 
intra- {U) and inter-orbital {U') Coulomb interaction as 
well as the Ising-component of the Hund's excha nge J 
{V = U—2J by symmetry; we set J = U/A as in [lolll^). 
The third line consists of the spin-flip contribution to the 
Hund's exchange (yielding together with the second line 
an SU(2)-symmetric contribution which can also be writ- 
ten as JSiiSi2, where Sim denotes the spin for orbital m 
and site i). The last term represents a pair-hopping term 
of same strength J. 

Method. QMC calculations which take the spin- 
flip component of Hund's exchange term into account 
have been a challenge. Although a straight-forward 
Hubbard- Stratonovich decoupling, exp{JATc\c2C^^C4) — 
(1/2) X]s=±i exp[s-\/ JAr(cJc2 — C3C4)], is possible, it has 
been recognized that it leads to a serious sign prob- 
lem Therefore, it was neglected in almost every 
DMFT(QMC) calculation so far, including 

To overcome this problem, several attempts have been 
made 0,0,0]. Among these, Sakai, RA, and Aoki 
proposed a new discrete transformation for the spin -flip 
contribution of the exchange and pair-hopping term[l^: 
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Here, A = i log(e2'^^^ + Ve*-^^^ - 1) , a = - log(cosh(A)), 
6EE;log(cosh(JAr)), fa = c\„C2a+cl„Cia, N^=ni„+n2a- 
2ni„n2a- The advantage of this decoupling is that the 
auxiliary field r is real in contrast to that of Hence, 
it is expected to yield better statistics in general [isj . 

However, even with this decoupling, we note that the 
usual Hirsch-Fye QMC algorithm does not work very 
well for DMFT calculation in the strong coupling regime 
or at low temperatures. For instance, for Hamiltonian 
and J = U /4, we found it to be infeasible to obtain 
a self-consistent DMFT solution for U > 2.2 when /3(= 
1/T) > 50 because the Green function G(t) has a large 
statistical error at r ~ (3/2. Therefore, it is difficult 
to clarify without ambiguity whether an orbital selective 
Mott transition indeed occurs in multi-orbital systems at 
low T by means of finite-temperature Hirsch-Fye QMC 
calculations; also see [iflj . 

Another recent advancement was the development of a 
new projective QMC (PQMC) algorithm by Feldbacher, 
KH, and Assaad In this algorithm, ground state 

expectation values (^'o|C'|^'o)/(^o|^o) of an arbitrary 



operator O are calculated as: 
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where Ht is an auxiliary Hamiltonian (its ground state 
\^t) is the trial wave function which is assumed to be 
non-orthogonal to the ground state |5'o) of H 14\). 

For Ht, it is convenient to take the one-body part 
of the Hamiltonian, because the limit (3 ^ 00 can be 
taken analytically in this case. Then, the starting point 
is the zero-temperature non-interacting Green function 
Gq{t, t') truncated to < r, r' < and discretized as an 
L X L matrix {L = 9/ At). From this Go{t, r'), the zero- 
temperature interacting Green function G{t,t') is ob- 
tained via the same updating equations for the auxiliary 
Hubbard- Stratonovich fields as for the finite-temperature 
Hirsch-Fye algorithm. 

While PQMC gives G(t) only for a limited number 
of (not too large) r-points, we need G{iuj) to solve the 
DMFT self-consistency cycle. To this end, the maximum 
entropy method (MEM) is employed to yield the spec- 
tral function A{u!) which allows for calculating G{iujn) — 

I '^'^ i^^^uj frequency iuJn- This makes a crucial 

difference to finite-temperature calculations. The large 
statistical errors occurring at t ~ /3/2 for finite temper- 
atures now occur for rather large r's. But even if there 
is a large statistical error for larger t's, the maximum 
entropy method can extract sufficient information from 
the first T points, discarding the larger t's with excessive 
statistical error. 

One of the main advantages of the PQMC method is 
that the convergence w.r.t. 9 is much faster than that 
w.r.t. (3 in the Hirsch-Fye algorithm 0]. (Note that the 
calculation time increases cubically for 9 and f3.) Hence, 
we take in the following PQMC calculations a finite 9 = 
20 (L = 64), which should be sufficiently close to the T = 
result: quantitatively small deviations are expected for 
larger 9's; qualitatively the behavior should not change 
anymore as in 0| . Similarly as in 13] , the central £ = 20 
are for measurement and V = {L — C)/2 ^22 time slices 
on the right and left side of the measuring interval for 
projection. Typically, we performed 7 x 10^ to 3 x 10^ 
QMC sweeps. 

Results. An indicator for the Mott-Hubbard transi- 
tion is the quasiparticle weight Z which is plotted in Fig. 
n^a). We clearly see that for the narrow band Z = 
for U > 2.6, while Z is still finite for the wide band. 
This means that there is a first Mott-Hubbard transi- 
tion in which only the narrow band becomes insulating 
at J7 ~ 2.5. This is consistent with the result of the 
DMFT(ED) calculation of [13, in which the critical Uc 
is estimated to be about 2.6. In contrast, there is a sin- 
gle first-order Mott-Hubbard transition at a smaller value 
[/c~2.1 if only the Ising-component of Hund's exchange 
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FIG. 1: (Color online) (a) Quasiparticle weight Z and (b) 
double occupancy D as a function of ?7 {J = U /4). Red (blue) 
squares (circles) are the data for the narrow (wide) band. For 
U = 2.4, two solutions are found: the wide band is metallic 
for both solutions whereas the narrow band is metallic (closed 
symbols) or insulating (open symbols). The solid triangle in 
(a) and (b) is the Uc estimate from DMFT(ED) [l^l; the inset 
enlarges the behavior around the transition. 

is taken into account (at r = 0.03; between J7c~l-8 and 
2.1 there are two coexisting solutions /hysteresis) jl^ . In 
our DMFT(PQMC) results, the wide band is still metal- 
lic at U = 2.7. But eventually, also the wide band has to 
become insulating at larger Coulomb interactions, since 
in the atomic limit both bands are insulating. (The cal- 
culation for larger Coulomb interactions unfortunately 
became computationally too expensive as even in the 
PQMC the statistical error brought about by the spin 
flip term of Hund's exchange increases dramatically.) 
Nonetheless, we can conclude from the data available 
that there are two different Mott-Hubbard transitions in 
which first the narrow and then the wide band become 
insulating. We have an orbital-selective Mott-Hubbard 
transition. 

In Fig.^b), the double occupancy D = (n^ni) for the 
two different bands is plotted as a function of U. We see 
that for the narrow band D « 0.01 for U > 2.6. A sim- 
ilar value of w 0.01 was reported Q for the one-band 
Hubbard model above the Mott-Hubbard transition, i.e., 
for U/W > 5.9/4. This suggests a Mott-Hubbard transi- 
tion very similar to the one-band Hubbard model, albeit 
only for the narrow band. 

Final evidence for the orbital-selective Mott-Hubbard 
transition is obtained from the spectral functions shown 
in Fig. 121 We can unambiguously say that the wide band 
is still metallic at [/ = 2.6, whereas the narrow band is al- 
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FIG. 2: (Color online) Spectral functions A{uj) for (a) the 
wide band and (b) the narrow band. For U = 2.6, the narrow 
band is insulating while the wide band is metallic. 

ready insulating with a pronounced gap. While the wide 
band shows a pseudo gap for an Ising-type of Hund's ex- 
change jl2j, our SU(2) symmetric result reveals a metallic 
peak in Fig. |21 

Let us now study the possibility of first-order Mott- 
Hubbard transitions. The first question is whether at 
U — 2.6 (where we find a metallic wide and insulating nar- 
row band) a second solution in which both bands are in- 
sulating (co)exists. Starting the DMFT self-consistency 
cycle with an insulating self-energy for the wide-band, we 
obtain however the very same (single) solution as in Figs, 
n and [3 This demonstrates that the orbital-selective 
Mott-Hubbard transition is not merged into a single first- 
order transition. There are two distinct Mott-Hubbard 
transitions. 

The second question is, Are the orbital-selective Mott- 
Hubbard transitions (generally) of first-order? In this 
case, two solutions should coexist for U < 2.6: one with a 
metallic and one with an insulating narrow band. Special 
care is necessary for insulating solutions in the PQMC 
with a very narrow charge gap. For such small charge 
gaps, might not be sufficiently large to project -via 
g-BH/2_ fj.Qj-^ ^YiQ metallic trial wave function onto the in- 
sulating ground state solution. We then note systematic 
errors even for intermediate r's, and substantial noise ap- 
pears in the charge gap of the spectrum calculated with 
the maximum entropy method. This makes the stabi- 
lization of a small-gap insulating solution delicate. This 
problem can be mitigated however by doing the max- 
imum entropy calculation with a reduced number of r 
points. Therefore, we used r points up to Tc ~ 3.2 and 
~ 1.6 for the following results. 

For almost all values of U, only a metallic or only an in- 
sulating solution is obtained for both rc~ 3.2 and ^1.6. 
However, for U = 2.4, we find both a solution with a 
metallic and with an insulating narrow band (the wide 
band is metallic in both solutions with only minor dif- 
ferences) . In Fig. 13 the spectral function of these two 
solutions are shown; the value of Z and D for the in- 
sulating solution is plotted in Fig. ^ as open circles and 
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FIG. 3: (Color online) Spectral functions A{uj) for (a) the 
wide and (b) the narrow band at U — 2.4. Two solutions with 
insulating/metallic narrow band coexist. 

squares. The DMFT(PQMC) data suggest that two so- 
lutions with metaUic and insulating narrow band coex- 
ists for U ^ 2.4, so that the Mott-Hubbard transition in 
which the narrow band becomes insulating (and in which 
the wide band stays metallic) is in general of first-order. 
Possibly, the insulating solution is metastable at r = 0. 

Discussion. For understanding the DMFT results it is 
instructive to remind ourselves of what is known for the 
two-orbital Anderson impurity model (AIM). If J > Tk 
(the AIM Kondo temperature) the impurity spins of the 
two orbitals form a steadfast spin-1 (triplet). For such an 
AIM and inequivalent orbitals it is known that this spin- 
1 is screened in two stages: first only by one orbital to a 
spin-i at T^, and then by the second orbital to a spin-0 
at r|. Within DMFT we now have to solve AIMs 

self-consistently: The AIM's T^'s of one DMFT iteration 
(crudely Tr- « ZW) sets the hybridization strength for 
the next DMFT iteration. Hence, we can interpret our 
DMFT results as following: Given the two inequivalent 
Kondo scales of the AIM, there is a t/-intcrval where only 
the hybridization strength (and Tk) of the narrow orbital 
is renormalized to zero by the DMFT iterations. Only 
the narrow band is insulating. 

If only the Ising-component of Hund's exchange is 
taken into account, the behavior of the AIM is completely 
different. Instead of a triplet, the impurity spins allign to 
5^ = ±1 (no Sz — component) . For J > Tk ( J ~ 0.5 and 
Tk^ZW^hOA at the Mott-Hubbard transition of [l^). 
there is no spin Kondo effect any more since it requires a 
spin-flip of the conduction electrons and, hence, a change 
of Sz by ±1. As soon as one orbital becomes insulating, 
there is also no orbital Kondo effect anymore: the whole 
system is unscreened, i.e., insulating. It is certainly in- 
teresting to study whether this kind of physics is relevant 
for magnetically anisotropic materials. 

Conclusion. Taking the full SU(2)-symmetry of 
Hund's exchange into account in the PQMC calcula- 
tion, we conclude that there are two consecutive Mott- 
Hubbard transitions, whereby -at least around the first 
transition- two solutions coexist. By clarifying the the- 



oretical side, we hope to stimulate further experiments 
on the orbital-selective Mott-Hubbard transition, e.g., in 
Sr^Ru04 where results were so-far negative in this re- 
spect 113 . 
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During the completion of our work, we learned about 
several related studies 0, . 
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